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ABSTRACT 


The Finite Element - Transfer Matrix (FETM) method has been developed to reduce 
the computations involved in analysis of structures. This widely accepted method, 
however, has certain limitations, and does not address the issues of control design. To 
overcome these, a modification of FETM method has been developed. The new method 
easily produces reduced models that are tailored toward subsequent control design. Other 
features of this method are its ability to (i) extract open loop frequencies and mode shapes 
with less computations, (ii) overcome limitations of the original FETM method, and (iii) 
simplify the design procedures for output feedback, constrained compensation, and 
decentralized control. 

This report presents the development of the new method, generation of reduced 
models by this method, their properties, and the role of these reduced models in control 
design. Examples are included to illustrate the methodology. 
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I. INTRODUCTION 


The straightforward applications of the existing theories to control of flexible 
structures is hindered by the dimension of their mathematical models. The control engineers 
have attempted to resolve this problem by resorting to some model reduction scheme. The 
control system that is designed based upon this reduced model is known to cause the 
detrimental "spill over" problems due to model errors, in general* . On the general topic of 
model reduction, many schemes have been developed which attempt to minimize model 
errors of one form or another. In the field of flexible structures, in addition to the structural 
analysts’ approach [2] whereby lower natural frequencies are retained, there is the truncation 
scheme based on performance sensitivity measures of the closed loop modes [3,4], modal 
cost analysis [5], and the residual flexibility procedure [6] which retains the static 
contributions of the deleted modes. With the exception of a few [3-5], all these schemes are 
based upon an invalid assumption that the modeling/model-reduction and the control design 
problems are separable. Once an attempt is made to integrate these two inseparable problems, 
it becomes evident that a good design philosophy can be generated only by employing the 
modeling techniques of structural analysts and the control theory of control engineers 
simultaneously. The need for control/structure interaction has been recognized as evidenced 
by special sessions in related conferences [7,8], and the first NASA/DOD Conference totally 
devoted to this subject [9]. Number of papers, in addition to Ref.[10], supporting this 
subject are contained in Ref.[ll]. Such an integration of the two disciplines forms the basic 
objective of this research effort. 

There exist several methods in computational structural mechanics that aid and 
improve the analysis and design of structures; Noor and Alturi offer a comprehensive survey 
of such methods in Ref.[12]. The methods that are of interest to the current investigation are 
those which produce lower order models; control engineers need these lower order models in 
order to overcome the "dimension problem" mentioned above. A method that seems to be 
well suited for this task is the finite element transfer matrix (FETM) method [13-20]. The 
FETM method has been developed to overcome the problem of large size matrices that are 
encountered while modeling flexible structures within acceptable accuracy. The main 
advantage of the FETM method is that it yields a reduced set of equations by operating at the 
(finite-) elemental level. The advantages offered by the FETM method are not fully utilized in 
the area of control design. This is since the FETM method has been developed primarily for 
the open loop analysis by structural analysts and has not yet reached the control community. 

* It has been pointed out in Ref.[l] that this may not always be the case. 
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The current research effort is an attempt to advance the FETM method, along with 
necessary modifications, to the control community as a viable tool for integrating structural 
modeling and control design phases. With a brief review of the FETM method and its 
limitations, Chapter II presents a new method that overcomes the limitations of the FETM 
method. Chapter n deals only with the extraction of natural frequencies and mode shapes. 
Chapter HI discusses the possible role of this new method in control design with some 
illustrative examples. The report is concluded in Chapter IV with recommendations and 
directions for future research. 
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II. A NEW METHOD BASED ON FETM 

The finite element technique is the most versatile tool in analyzing dynamic response 
problems of complex structures, particularly in evaluating natural frequencies and mode 
shapes. In order to describe a structure with good accuracy, it is necessary to divide the 
structure into a very large number of elements which in turn necessitates a large number of 
degrees of freedom (dof). This leads to large mass and stiffness matrices requiring significant 
amount of computer time to process. Consequently, alternative methods for reduction of the 
number of dof associated with a given structure have been developed. One of the most effective 
methods of order reduction which appears to be compatible with active control is the finite 
element transfer matrix (FETM) method [13], 

The transfer matrix method is based on the selection of a state vector which is 
transmitted across the structure upon multiplication by a transfer matrix - the state vector 
consists of the dof and the generalized nodal forces along a section. Realizing that the number 
of dof in a finite element analysis of a structure can be considerably reduced by combining with 
the transfer matrix approach (which is most efficient when it is used along a dominant 
direction of a 'chain-like' structure), Dokainish [13] was the first to develop a combined 
method which exploited the best features of both the methods. He was able to show that the 
FETM method led to a frequency determinant of size 18 by 18, while the standard finite 
element technique produced a 108 by 108 determinant, the number of nodes being the same in 
both cases. The main advantage of the FETM method is that it yields a reduced set of equations 
by operating at the elemental level. 

There are two limitations of the FETM method of Dokainish : (1) the method is 
applicable only to 'chain like’ structures (these include beams and rectangular plates, but not 
circular plates), and (2) it requires the number of dof along the edges of all substructures to be 
the same. Several authors [14-20] have proposed means to overcome the second limitation by 
incorporating certain approximations. In these modified methods, though the second limitation 
of the FETM method is reduced, it is not completely eliminated. In what follows, we briefly 
discuss the FETM method of Dokainish, and present our new method, and conclude this 
chapter with an example. 


2.1 Brief Review of FETM method 

The transfer matrix method is generally associated with the concept of stiffness. 
Pestel and Lechie [21] demonstrated that for one-dimensional structural members such as 
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beams, the transfer matrix could be derived from the stiffness matrix. Dokainish [13] was 
the first to suggest a combined FETM method. He recognized that the transfer matrix method 
could be extended to two-dimensional chain like structures by deriving the transfer matrix for 
a plate strip from the corresponding stiffness matrix. Since the publication of Dokainish's 
paper, several authors (notably [14-20]) have proposed refinements and extensions of the 
FETM method. All of these authors derive the transfer matrix from the stiffness matrix in 
much the same fashion as Dokainish. This derivation, following that of Ref. [20], is 
summarized below. 

The equation of motion of a structure can be written using a dynamic stiffness matrix 

S e as 

S A = f; S e = K e - ; ■ (2.1) 

where Ke is the stiffness matrix. Me is the mass matrix, and to is a natural frequency of the 
structure. Here q e are the nodal displacements and f are the nodal forces. The standard finite 
element technique extracts the natural frequency u> by setting ISI — 0, where S is a submatrix 
of S e , obtained after enforcing the boundary conditions that require some of the dof to be 
zero. Note that the size of this determinant is as large as the number of the total free dof used 
to describe the structure. 

Consider now a structure which consists of N substructures. The dynamic stiffness 
equations of each structure can be obtained as in Eqn.(2.1) by employing finite element 
techniques. These equations for the ith substructure can be partitioned as 

s; sLlii 

S H Kli 

where f denotes the internal generalized forces; the subscripts T and V refer to quantities on 
the left and right boundaries of the substructure. Note the single subscripts on the diagonal 
submatrices Sj 1 and S^. In writing Eqn.(2.2) we have assumed that either there are no interior 

nodes or the interior nodal displacements have been eliminated as mentioned in Ref.[21]. Hie 
transfer matrix T* for the ith substructure is then found by manipulating Eqn.(2.2) to yield 

-isiA* [sii 1 (<■; 

.-s^is'r-s; -sjsLr'lK 

or simply 


-t 


(2.3a) 


( 2 . 2 ) 
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A = *; • 


(2.3b) 


Note the definition of the state vector xj at the right boundary. Now for the adjacent (i+l)th 
substructure, equilibrium and continuity of displacements require that 

(2.4) 


i+l i 
X, = x r 


Equations (2.3b) and (2.4) are then combined to give 
i+l ™4+l r^i i 


x ; +l = rfij . 


(2.5) 


Extending Eqn.(2.5) to N substructures results in 

x" = TxJ ; T = T^T*' 1 ... T 1 , (2.6) 

where x^ and Xj 1 are the state vectors on the extreme right and left boundaries of the entire 
structure. Enforcing appropriate boundary conditions on x^ and x{ would require ITfl = 0 


for a nontrivial response, where Tf is a submatrix of T. This condition is then used to extract 
the natural frequencies of the structure. Note that the size of this determinant is equal to the 
number of dof along a boundary. 

The derivation of the transfer matrix from the dynamic stiffness matrix requires the 
inversion of the submatrix as shown in Eqn.(2.3a). In a strict sense, this inversion is 

possible if and only if S^. is a square nonsingular matrix. However, Sj, is square if and only 


if there is an equal number of nodes (actual dof) on the right and left boundaries of the 
substructure, which in general may not hold. The submatrix , thus, cannot be inverted in 

general. Therefore, the above formulation of combined FETM method is only applicable to 
models which have the same number of dof on all substructure boundaries. 

Recognizing that rectangular transfer matrices could occur in practice, Pestel [19] 
proposes using a left-inverse of . This is possible only if the number of nodes on the left 

boundary of a substructure is equal to or greater than that on the right boundary. This 
restricts Pestel's approach to structural models in which the number of nodes on the 
substructure boundaries decreases (or increases) monotonically from one exterior boundary 
to other. As suggested by Degen, et al. [20], a further extension of Pestel's procedure to 
rectangular transfer matrices would be to use a generalized-inverse, rather than a left-inverse, 
of S},. . This however introduces approximations into the derivation of the transfer matrices 

from the stiffness matrices. More recent developments to overcome these limitations have 
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been reported be Degen, et al. [20], employing a mixed finite element formulation. In [20], 
the energy expression is obtained as Reissner functional, which allows the governing 
equations to be obtained in mixed form, that is, as a combination of stresses and 
displacements. This formulation permits the boundaries of a substructure to have a different 
number of dof "within certain limits"; the limits being that the number of nodes on the right 
boundary be not "too large" relative to the number on the left boundary. Hence, the limitation 
of the original FETM method, though reduced, are not completely eliminated. 


2.2 A New Method 


In this section, we present a new method, which completely eliminates the above 
limitations, without introducing any approximations. The development of this method is 
based on the fact that the natural frequency of a structure must render ISI = 0, where S is the 
dynamic stiffness matrix of the total structural model that results after enforcing the boundary 
conditions. That is, for an arbitrary partitioning of S, the natural frequency must satisfy 


S 


l 




(2.7a) 


which, under the assumption of nonsingular Si yields 

|s 2 -s 21 s; l s 12 | = o. 


(2.7b) 


In what follows, we obtain Eqn.(2.7b) from operating at substructure level, instead of 
beginning with Eqn.(2.7a) which in general is of large dimension. 


2.2.1 Development 


Consider Eqn.(2.2), which is rewritten below for convenience 

,i,T 


s i 

S L 

» - 

i 




s‘ 

r 

i 

[%\ 


< 


; s; = [S 1 *]' 


( 2 . 8 ) 


where q represents the 'free' dof (i.e., after enforcing the boundary conditions) on the left 
and right edges of the substructure as indicated by the subscripts T and V respectively. 
Using the top part of Eqn.(2.8), 


6 



(2.9) 


F U = I 

r U A ’ 


the qj dof can be eliminated to obtain a reduced representation of the ith substructure in terms 
of only qj: 

s r.i<u = (210a) 


where. 


c i.i _ rc* , H^' — Q* o* rQ*'V^Q* • c’’’ — c* 

S f i - [S f J \ ^ » ^1,1 “ *i * 

Tji>i A rcWi'iij'i 

F r,l = ' S rlt S l,lJ F U 


(2.10b) 

(2.10c) 


Note that the inverse involved in Eqns.(2.10b,c) is only of a square symmetric matrix; we 
assume the existence of the inverse* . The V in the subscripts and superscripts are used to 
indicate that the corresponding elements ( S and F in Eqns.(2.10a,b,c)) include some 
’transferred’ information - the information transferred is indicated by the subscripts and 
superscripts following the 

Now consider the dynamic stiffness matrix of the (i+l)th substructure: 


(2.11) 


Under the assumption that there are no external forces present along the edge common to the 
ith and the (i+l)th substructures, using Eqn.(2.10a) and the equilibrium condition 

fi* l + f| = 0, (2.12) 


sr 

_i+l 

i+l' 

9l 


Y 

c i+l 
. rl 

sj+t 

i+1 

% 


Y. 


renders the top part of Eqn.(2.1 1) as 

_i+l i+1 , c i+l i+1 t-M A c i.i i 

s , q, +s i, % - • 

In view of the continuity condition 

i+i i 
9i “ 9 r » 

Eqn.(2.13) becomes 


j+l.i 


c i+l,i i+i Q i+1 i+l _ t?+ 

S l,l 9 l + S lr - F l.l 

where, using the notation introduced in Eqns.(2.10b,c), 


(2.13) 


(2.14) 


(2.15a) 


* The significance of this assumption will be discussed later in Section 2.2.2. 
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s$ w -0 T 4sr*s“, 


(2.15b) 


pi + l.i A pU 
1.1 “ r.l ‘ 


(2.15c) 


Note that, Eqn.(2.15a) is similar to Eqn.(2.9). Therefore, one can proceed as before to 
eliminate q, 1+ from Eqn.(2.11). This yields a reduced representation of the ith and (i+l)th 

substructures in terms of only q£ +1 : 

c i+l,i i+l 4+1 -j+l, i 4 

S ,4 - r + F ,J V (2.16a) 

where, 

c* +1>i _ A «i+l oi+l-ni+l.L-l-i+l 

\i ~ LS r,i J - S r - S rl [S u ] S h , (2.16b) 

r?i+l,i A c i+l re i+l,i--l T a+l,i 

F rj * - S rt t S l,l I +1,1 • ( 2 . 16 c) 


Proceeding similarly, a reduced representation of the substructures i, i+l, i+2, .... j 
(j ^ i), in terms of only qj can be obtained as 

• (2.17a) 

where, for k = i,i+l, ... j-lj, 

s“ = [sft T A s t k .sXj]"s^ 

S U - [^) T - S k + S k , u ; s|j’‘ A 0 


(2.17b) 


pa., . i £ u 5 k >i 

; k = i 


{: 


We emphasize that the matrix inversions needed in Eqn.(2.17b) are only of square 
matrices. Hence, as long as these matrices are nonsingular, Eqn.(2.17a) can be obtained 
without introducing any approximations. Moreover, no assumptions have been made 
regarding the number of dof on the left and right edges of the substructures - the boundary 
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conditions of the substructures can be arbitrary. Thus, the second limitation of Dokainish's 
FETM method is totally eliminated. Means of overcoming the first limitation is described in 
Section 2.2.3. 

2.2.2 Algorithm 

We illustrate the use of Eqn.(2.17) in determining a natural frequency (co) of a 
structure that is divided into N substructures - the substructures are numbered from left to 
right It is assumed that the mass and stiffness matrices are available for all the substructures, 
hence, the elemental dynamic stiffness matrices can be calculated. We propose the following 
algorithm. 


Determination of natural frequency 


StepO. 
Step 1. 


Assume an initial value for co. 
Set 



k = 1 


Step 2. 


a) Calculate 

Q k,l _ c k Q k-l.l 
\l - a l + ^r.1 


(2.18a) 


Step 3. 


r.l “ r ‘ rlL\l J b lr 

b) If k = N, go to Step 3. 

If k < N, set k = k + 1, and repeat Step 2. 


a) Determine 

A<*» - |tf| • 


(2.18b) 


(2.19) 


b) If IA(co)l £ e, a prespecified representation of numerical zero, stop. The value of 
co is a natural frequency of the structure. 

If not, use an "appropriate" value for co and go to Step 1. 
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The justification of this algorithm is based on the following theorem. 

Theorem 2.1 

Under the assumption of nonsingular for 1 £ k £ N, the natural frequency (O of the 
structure satisfies A(co) = 0. 


Proof: 


Let N = 1. Then CD satisfies 


S 1 S 1 
5 1 nr 


s 1 s 1 

rl 5 r 


- MM - o 


( 2 . 20 ) 


Since S* = S^j 1 is assumed to be nonsingular, Eqn.(2.20) results in, for N = 1, 

Ato) = Isfj'l » o. 


(2.21) 


Now let N = 2. Then co must satisfy 


0 


s i sVsf s l 


0 K K 


= Is 


S 2 + S lfl 


- | s .1| s u || s ?.'| = 0. 


which again yields Eqn.(2.21) for N = 2, under the assumption of nonsingular Sj 1 and SjY 
Proceeding similarly for a general N proves the theorem. # 


2.2.3 Remarks on the Algorithm 


A few remarks pertaining to this algorithm are in order. 
a) The nonsingularity assumption ofsff 

For the above algorithm to work it is necessary for sjj 1 (1 £ k £ N) to be 

nonsingular - an assumption made in the development of the algorithm. The significance of 
this assumption follows. Define 

8 k «0) = |sf,j . (2.22) 

Then S^ 1 is singular if and only if 8k(co) = 0. Let k = 1. Then 
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(2.23) 


SjCco) = 



where S* is the dynamic stiffness matrix of the first substructure with its right edge fixed, 
i.e., qj = 0. Therefore, 8 i(cd) = 0 implies that the assumed frequency is a natural 


frequency of the first substructure with its right edge fixed. Hence, the nonsingularity 
assumption of S^J 1 only implies that the assumed value of © is not one of the frequencies of 


the first substructure with its right edge fixed. We shall assume this is the case, and let 
k = 2. Then 


5 2 (o» = |S?*| = l^l-lsX-silS^si 


V©) 


i.i 


U 

i 


sVsf 


(2.24a) 


Consider now a structure that consists of only the first and second substructures with the 

2 

right end of the second substructure being fixed, i.e., q,. = 0. Then the natural frequency of 


this structure must satisfy 


sWsf 


- 8j(©) 8 2 (©> = 0. 


(2.24b) 


Therefore, comparison of Eqns.(2.25a) and (2.25b) reveals that f> 2 (©) = 0 if © is a natural 
frequency of this "2-element" structure whose right edge is fixed. Proceeding similarly, we 
arrive at the conclusion that the nonsingularity assumption of S£j is violated if the assumed 


value of © is one of the frequencies of a structure that consists of only k substructures with 
its right edge fixed. 


b) The choice of an "appropriate" (0 

Note that in Step. 3b there is a need for an appropriate updating of ©. Currently, we 
are using a variation of the successive bisection method to update ©, and it seems to work 
efficiently for the beam examples considered herein. There may exist other methods that 
work more efficiendy . An appropriate method needed in Step 3b is under investigation. 
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c) Inclusion of spring-support/ point-mass 

The presence of any spring-support and/or concentrated point masses at any of the 
intermediate nodes of the structure can be easily accommodated in this algorithm. This can be 
achieved as follows. Suppose that some nodes along the edge common to ith and (i+l)th 
substructures have spring supports and/or point masses. Then, to take these into account, the 
equilibrium condition Eqn.(2.12) must be expressed as 

fj +1 + = IK!,, - ®V W ] qi" (2.25a) 


where Ki (resp. Mj +1 ) is a diagonal matrix consisting of the 




supports (resp. masses of the point masses) • some of these may be zeros along its diagonal. 
Then, the dynamic stiffness matrix J^’Wuld correspondingly change to 


S U 1J - 


(2.25b) 


The presence of spring-support and/or point-masses along the left and right boundaries of the 
overall structure can be similarly taken into account 


d) Determination of Mode Shapes 


In the absence of any external forces, Eqn.(2.17a), for j = N and i = 1 yields 

c N,i „ N o 
S r,l % m °» 


(2.26) 


from which, in view of A(co) = 0, a nontrivial q^ can be determined. Then using the 

equation, analogous to Eqn.(2.15a), 
oN. 1 N.« N N n 

S U =° < 2 - 27 ) 

N N 1 

qj ( = q r ' ) can be determined. Proceeding similarly, the entire mode shape of the structure, 
corresponding to that frequency, can be determined. 


e) Development from assembled equations. 

Note that the computations involved in Step2a are repetitive, i.e., they are repeated 
for k = 1,2, ... ,N, indicating that the computations are carried out at substructure level as 
one proceeds from the first (k=l) to the last (k=N) substructure via the intermediate 
substructures (k=2,3, ... ,N-1). The final outcome of this step is the matrix S^’ 1 . This matrix 
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can also be obtained from the complete assembled dynamic stiffness matrix S representing 
the total structure, though this involves inversion of larger matrices. This follows from the 
proof of Theorem2. 1 . 

2.2.4 Variations of the Algorithm 

Variation A. 

Note that the development, and the algorithm, assume that in obtaining Eqn.(2.19) 
one advances from left-to-right as the computations are performed at the substructures level. 
One could also advance from right-to-left to determine the natural frequency. 

Determination of natural frequency: Variation A 

Step 0. 

Step 1. 

Step 2. 


Step 3. 


This approach amounts to transferring all dof (c£; k = N,N-1, .. . ,2) to q^ in order to obtain, 
as in Eqn.(2.26), a reduced representation of the structure in terms of qj 1 : 


Assume an initial value for ©. 
Set 

sf; 1>N = 0; k = N 


a) Calculate 

c k.N _ c k k+ 1.N 
^r.r ~ b t + b U 

pk|N nk rtk|*Qk,N -.-l^k 

^l,r " * V\r J b tl 

b) If k = 1, go to Step 3. 

If k > 1, set k = k - 1, and repeat Step 2. 


(2.28a) 

(2.28b) 


a) Determine 
A(0>) £ Is;/ | 


(2.29) 


b) If IA(co)l £ e, a prespecified representation of numerical zero, stop. The value of 
© is a natural frequency of the structure. 

If not, use an "appropriate" value for © and go to Step 1. 
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C 1.N 


1 

<k 


= 0. 


(2.30) 


In fact, one can obtain such a representation of the structure in terms of any dof, qj ( or qj ) 


as follows: 


Variation B. 


Proceed from left-to-right to obtain 



and proceed from right-to-left to obtain 


C i+1,N 

b l.r 




Now, enforcing the equilibrium and continuity conditions yields 

re ** 1 i ni+l|N. i+l a 
t S r,l +S l,r H =° 


The Eqn.(2.19) will then be replaced by 


A(co) = 


s 1 ;! + s* +1,N 

r,l lr 


(2.31a) 

(2.31b) 

(2.31c) 

(2.3 Id) 


This procedure is summarized below: 


Determination of natural frequency: Variation B 

Step 0. Assume an initial value for to. 

Phase I. 


Step 1. 1. Set 


, 0,1 


Step 1.2. 


KS = °; k» i 


a) Calculate 

c k>1 _ c k , c k - l,l 

5 1,1 ~ b l + b r,l 

C*k,l ^ c k nk r okfl-i-l pk 

r,l ~ b r " rl^l.l J b lr 

b) If k = i, go to Step 3. 

If k < i, set k = k + 1, and repeat Step 1.2. 


(2.32a) 

(2.32b) 
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Ph ase I J, 

Step II. 1. Set 

S?* 1 *- 0; k = N 

Step II. 2. 

a) Calculate 

„k,N k Q k+1,N 
^r.r = + \r 


(2.33a) 


Step 3. 


k.N _ „k (jk .kN.-l k 

^l,r ” ^1 ' \^T,T J b Tl 

b) If k = i+1, go to Step 3. 

If k > i+1, set k = k - 1, and repeat Step 112. 


(2.33b) 


a) Determine 

a(co) = |s;j + s;; 1N | 


(2.34) 


b) If IA(cd)I <i e, a prespecified representation of numerical zero, stop. The value of 
co is a natural frequency of the structure. 

If not, use an "appropriate" value for co and go to Step 1. 


Variation C. 

Moreover, the algorithm fen* determining the natural frequency can also be obtained in 
terms of any combination of the dof. Suppose, for example, that we wish to develop this 
algorithm in terms and qj (1 < i < j < N). Then, we could proceed as follows. Proceed 


from first substructure to ith substructure to obtain 

S U <i + S U = °- (2.32) 

Proceed from Nth substructure to jth substructure to obtain 

^♦S-qj-O. (2.33) 

Now. to eliminate the intermediate dof, we begin with (for convenience) the dof common to 
ith and (i+l)th substructures: 
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(2.34a) 


s U; + t s > s i +1 K i+s L +1 q; +1 - o* 

and compute 

i rc M+l-,-l t ni i ot+1 i+1 . 

q r = -[S r ] { S rl q 1 + S ]r q f } , 


(2.34b) 


where 


„i,i+l A c i c i+l 

&r - + • 

(2.34c) 

Substituting Eqn.(2.34c) in Eqn.(2.32) yields 


[S ]4 * Sfr(S* ,i+1 ) 's*,] qj - S^(si’ i+1 )* 1 sj^’ 1 q* +1 = 0. 

(2.35) 

Now consider the equations for q r l+1 from the (i+l)th substructure: 


c i+l i+1 c i+l i+1 -i+1 

S rl <*1 +S r % =f r * 

(2.36a) 

which in view of Eqn.(2.34b) and since qj +1 = qj. , becomes 


-S r i (S r ) S rl qj + [S r -S ri (S r ) s t ]q f = fj . 

(2.36b) 


Then, using the equilibrium condition, q r l+1 can be expressed in terms of qj and qj +2 which 
expression can then be used to write Eqn.(2.35) in terms of qj and q‘ r +2 Proceeding 


similarly the following equations can be developed: 




■ ■ 
1 

*ii 

*12 


*21 

*22 

i 


(2.37) 


where Sij are obtained through the above sequential calculations. Eqn.(2.37) can then be 
used to determine the natural frequency. 

Even though we have offered the above development - from Eqn.(2.32) to 
Eqn.(2.37) - as a variation of the algorithm of Section 2.2.3, this algorithm, along with its 
variations A and B are really special cases of Variation C. We have chosen to offer Variation 
C last in order to make the development more comprehensible. 

Observe that Variation C of the algorithm can be used for structures that are not chain- 
like, such as circular plates, thereby eliminating the second limitation of the Dokainish's 
method. 
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2.3 Example 


A simply supported beam, with some minor variations, is employed in this section to 
illustrate the applicability of the new method proposed herein. The beam under 
consideration is shown in Fig.2.1 below. 


1 A 2 B 3 C 4 



< 2IV3 »| 

◄ L * 


Fig. 2.1 Simply Supported Beam: 3 Elements 


The beam is modeled by three (hence, N = 3) cubic beam finite elements A,B, and C with 
points 1,2,3, and 4 as the nodal points. The linear and the angular deflections of the beam 
are denoted by x and 0, respectively. The numerical value assumed for this beam are: the 
length of the beam L = 1; pA = 420, where p is the mass density per unit area, and A is 
the cross sectional area; El = 1, where E is the modulus of elasticity, and I is the bending 
moment of inertia. 

The following method illustrates the new method, as applied to this example. 


The elemental mass and stiffness matrices are given below, with the elements 
numbered from left to right 



‘ 1.4815e-l 

-1.4444e+0 

-l.lllle-1 


1.2000e+l 

5.4000e+l 

6.0000e+0 

M 1 = 

-1.4444e+0 

5.2000e+l 

2.4444e+0 

;k'= 

5.4000e+l 

3.2400e+2 

5.4000e+l 


-l.lllle-1 

2.4444e+0 

1.4815e-l 


6.0000e+0 

5.4000e+l 

1.2000e+l 



5.2000e+l 

-2.4444e+0 

1.8000e+l 

1.4444e+0 

-2.4444e+0 

1.4815e-l 

-1.4444e+0 

-l.lllle-1 

1.8000e+l 

-1.4444e+0 

5.2000e+l 

2.4444e+0 

1.4444e+0 

-l.lllle-1 

2.4444e+0 

1.4815e-l 
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3.2400e+2 

-5.4000e+l 

-3.2400e+2 

-5.4000e+l 


-5.4000e+l 

1.2000e+l 

5.4000e+l 

6.0000e+0 


-3.2400e+2 

5.4000e+l 

3.2400e+2 

5.4000e+l 


-5.4000e+l 

6.0000e+0 

5.4000e+l 

1.2000e+l 



' 5.2000e+l 

-2.4444e+0 

1.8000e+l ‘ 


' 3.2400e+2 

-5.4000e+l 

-3.2400e+2 

M*= 

-2.4444e+0 

1.4815e-l 

-1.4444e+0 

;K 3 = 

-5.4000e+l 

1.2000e+l 

5.4000e+l 


1.8000e+l 

-1.4444e+0 

5.2000e+l 


-3.200e+2 

5.4000e+l 

3.2400e+2 


Hence, for the first element, the dynamic stiffness matrix is obtained as 
sj = 1.2000e+l - 1.4815e-l * to 2 

S i = t S ril T =[ 5 -4000e+l+1.4444e+0*to 2 6.0000e+0+l.lllle-l*a) 2 ] 

3.2400e+2-5.2000e+l*o) 2 5.4000e+l-2.4444e+0*co 2 

s r = 

5 .4000e+ 1 -2.4444e+0*co 2 1.2000efl-1.4815e-l*o> 2 

Similarly the dynamic stiffness matrices for other elements can also be obtained. With these 
matrices available we are now in a position to use the algorithm of Section 2.2.2. 

Using 0.48184 as the initial value for co yields A(co) = 0.0047. The variation of A(o>) with 
to is shown in Fig.2.2. The Table 2.1 gives the frequencies obtained by conventional finite 
element method and by the new method, and the value of A(co) corresponding to these 
values of the frequencies by the new method. 


Table 2. 1 Results for the system of Fig.2. 1 


co by finite element 

to by new method 

A(to) 

0.48184 

0.48184 

0.0047 

1.94960 

1.94960 

0.0020 

4.81120 

4.81120 

0.0043 

8.94700 

8.94470 

0.0015 

16.0110 

16.0110 

0.0021 

22.0580 

22.0580 

0.0044 
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») .vs. co for the system of Fig.2.1 


In the following few pages, we present the results obtsined by employing the new 
method to different configurations of the simple beam. The numerical value used for the 
stiffness of the spring in Fig. 2.7 is 15 units. Note that the original FETM method can not 
applid to the configurations of Figures 2.5 and 2.7. 

1A 2 B 3 C 4 D 5 



co by finite element 

co by new method 

A(co) 

0.48380 

0.48380 

0.0043 

1.93430 

1.93430 

0.0070 

4.41200 

4.41250 

0.0278 

8.53690 

8.53740 

0.0214 

13.5670 

13.5900 

0.0270 

21.4970 

21.4970 

0.0004 

31.9440 

32.2200 

0.0028 

38.7870 

39.2270 

0.0040 



co by finite element | co by new method I A(co) 


0.54309 0.54280 0.0500 


0.87665 0.87600 0.0240 


2.22700 2.22560 I 3.7967 


3.13510 3.12350 I -0.6428 


5.70690 5.70870 0.0270 


8.25750 8.25876 0.0004 


10.3830 10.5223 J 0.0028 




























































1A 2 B 3C4D5 



Fig. 2.7 spring supported, simple beam 


Table 2.4 Results for the system of Fig.2.7 


co by finite element 

co by new method 

A(<o) 

0.17611 

0.22100 

0.0012 

0.47222 

0.48600 

0.0097 

1.13360 

1.12100 

0.0324 

2.13720 

2.13760 

0.0257 

3.40080 

3.40400 

0.0343 

5.37510 

5.37510 

0.0026 

8.05350 

8.05350 

0.0233 

9.79780 

9.79780 

0.0421 


23 































id J u ^ 



III. APPLICATION TO CONTROL DESIGN 


The FETM method [13], along with its modifications [14-20], and the new method 
presented in Chapter n, have been developed for open loop (free vibration) analysis. A 
control engineer, however, is more concerned with the closed loop design and analysis, and 
requires an appropriate reduced model of the structure. Models consisting of finite 
dimensional ordinary differential equations are preferable to those consisting of partial 
differential equations, on the basis of practical considerations, such as the implementation of 
the control algorithms. Traditionally, the finite element technique is adopted to produce such 
models. These models, for structures with neither dissipative nor gyroscopic elements, have 
the form 

Mq + Kq = Eu ; qeR n , ueR m (3.1) 

where the second-order state vector q consists of all the degrees of freedom, and the input 
vector u represents the control inputs (actuators) to the system. The symmetric mass matrix 
M is assumed to be positive definite, while the symmetric stiffness matrix K is assumed to be 
nonnegative definite. 

Since n, the dimension of q, equals the total degrees of freedom, it can be very large - 
the larger the number of substructures that the overall structure is divided into (as is 
necessary to render the model of Eqn.(3.1) to be a more accurate representation of the 
structure), the larger is n. As a consequence, the dimension of the model hinders the 
straightforward applications of the existing control theories. Therefore, a need to represent 
the given structure by lower order models of acceptable accuracy arises. 

In the field of flexible structures, the popular model reduction schemes are : the 
structural analysts approach [2] which retains modes corresponding to lower natural 
frequency; the modal cost analysis [5] which retains modes that have the largest contribution 
to the mission objective, measured via quadratic cost function; and the open- and closed- loop 
balancing method [22] which produces reduced models based on a "measure of importance" 
of the modes. 

Note that the above model reduction methods operate at ’modal’ level - i.e., they 
produce reduced models based on the modal data. Hence, the starting points for these 
methods rely upon the availability of reliable modal data. This is difficult to obtain for 
Eqn.(3.1), particularly when n is large - it is precisely due to this difficulty that the FETM 
method has been developed. 
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In this chapter, we show that the new method developed in Chapter II offers a 
procedure to aid in the control design phase. The advantages of this new method as applied to 
control design of flexible structures are : 

a. it does not require the modal data. It needs only the frequencies which , as shown in 
Chapter n, can be obtained in passe. 

b. it operates at the elemental/substructures level; hence, it does not require even the 
assembled equation (3.1). 

c. the reduced representation produced by this method can be tailored to suit the control 
architecture, such as sensor/actuator locations, centralized/decentralized feedback, 
output/state feedback, and mission objectives. 

The properties of the resulting controllers are yet to be investigated. 

3.1 Development 

In order to establish that the procedure developed in Chapter II can be easily adapted 
for control design, we begin by considering the homogeneous version of Eqn.(3.1): 


Mq+Kq = 0 ; qeR n . 

(3.2) 

As usual, assuming a harmonic response, we get 


GO 

§ 

II 

* 

S' 

A 

II 

O 

* 

(3.3a) 

which in partitioned form, with ni being arbitrary, is written as 


Sj qj(t) + S 12 q 2 (t) = 0, q t eR 1 

(3.3b) 

S T n qfi) + S 2 q 2 (t) =0, q 2 6 r" 2 ; n 2 =n-n x 

(3.3c) 

where 


Sj = Kj-cd 2 !^ , S 2 = Kj-G) 2 !^ , S 12 = K l2 -(0 2 M n . 

(3.3d) 

Using Eqn.(3.3b) we obtain 


q^t) = - S 12 q 2 (t) = " Sj S 12 q 2 (t) 

(3.4) 

establishing a relationship between qi(t) and q 2 (t). Hence, the lower part of Eqn.(3.2) yields 

M 2 1 q 2 + K 2j q 2 - 0, 

(3.5a) 

where 


^2,1 = ^*2 ’^12^1.2’ ^2,1 = ^2 “^12® 1,2 ' 

(3.5b) 
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Note that Eqn.(3.5a) is of dimension n 2 , and it still represents the original 'large' model (3.2) 
of dimension n with the dof q* "absorbed" into q 2 . Also note that the relation (3.4) and the 
model (3.5a) are accurate for a homogeneous system. If the system is not homogeneous, 
these equations are only approximate, with the degree of approximation depending on the 
forcing function. 

Equation (3.2) represents the assembled set of equations of motion describing all the 
substructures, and hence is of large order, that is, n may be quite large. Since n is large, ni 
will still be large, particularly when n 2 is required to be small. If nj is large, then the matix 
being inverted in Eqn.(3.4) is large, which may cause numerical difficulties. In order to 
avoid this difficulty, recall that this approach, as presented thus far, operates only on the 
assembled set of equations, and not at elemental level as developed in Chapter II. Therefore, 
we would like to develop a method that operates at elemental level - and hence, requires less 
computations - which would produce the same representation (3.5a). This development is 
given below. 

Consider an ith substructure and proceed analogous to the development in Section 
2.1, but involving mass and stiffness matrices in place of the dynamic stiffness matrix.This 
would yield, in terms of accelerations and positions, in place of Eqn.(2.9) 




where 


Ka 4 K- Muir's;, ; 




(3.6a) 


(3.6b) 


and where the superscripts and subscripts follow the same notation adopted in 
Eqns.(2. 10a,b). Note that the matrices and Kjlj 1 are no longer symmetric, in general. 


Now consider the (i+l)th substructure, and suppose that there are some external 
(generalized) forces u> applied along the interface of the ith and (i+l)th substructures. Then 
the equilibrium condition (2.12) becomes 

fj + fj +1 = u* . (3.7) 

(If there are any concentrated masses and/or spring supports that exist along the interface, 
they can also be accommodated as pointed out in Remark (c) of the previous chapter). Use 
Eqn.(3.7) and the continuity condition (2.14) to obtain the equations of motion for the 
interfacing dof as 
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r _ i+1 -i,L •• i+1 . rir i+l t ^ i,i, i+1 , *i+l •• i+1 , x^i+1 i+1 

[M; +1VTJ qj + [Kj +K 1 lq 1 +M^ q r + q. 


,4+1, i i ,4+1, i+1 i 

= *1,1 h + *ti u 


(3.8) 


Continuing again with the assumption of harmonic response, the equations of motion 


of the ith and (i+l)th substructures assembled together, can now be written in terms of q,. 
as 

i+l,i -i+1 -4+14 i+1 4+1 -4+14 4 -4+14+1 


i+1 


» -1+1,1 -1+1 --1+1,1 1+1 ,4+1 pl+1,1 J -4+1.1+1 1 

K,1 % +K r.l <4 = f r + F r,l + *V,1 U 


(3.9a) 


where 


_ -i+l,i A . -i+1 » -i+1 |.-,i+l,i-.-l ^ i+1 „i+l,i A «i+l „i+l--,i+l t i--l-,i+l 

r,l =M r " M rl t S U 1 S lr ’ K r,l = K r “ K rl t S U 1 S lr 


(3.9b) 


_i+l,i+l A t ,i+l re i+l,i+l 1 -l_i+l,i+l 
F r,l “ S il t S l,l J F l,l 


Thus, the development presented thus far offers a procedure to obtain a reduced 
representation of the structure consisting of the ith and (i+l)th substructures. If, for example, 
the total structure consists of only these two substructures (i.e., N=2, and i=l), and if there 
are applied forces u° along the left boundary (making f*=u°), u 1 along the interface, and u 2 

along the right boundary (making f^=u 2 ), then from Eqn.(3.9a), a reduced model of this 


structure would be 

X* 2 , 1 "2 „2,1 2 ^,1 2,1 

Kj <1, + K ,,i <i, - <1 u ,j • 

where 

[<‘] T = A [(uY,(u') T . ( uy] T 
= t ^ . ^, 2 . 1 1 


(3.10a) 


(3.10b) 


These equations can be derived directly from the results of Chapter n. This statement is 
verified by recognizing that 

which provides the "mass” and "stiffness" matrices required of the reduced models. 
Moreover, from the discussions on Section 2.2.4, Variations of the Algorithm, such a 
reduced model can be developed in terms of any selected dof - i.e., in terms of dof either on 
the left boundary, or on the right boundary, or on any intermediate interfaces, or any 
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combination. This flexibility is required in obtaining models that are to be used for 
subsequent control design, since the dof we wish to keep in the reduced models will differ in 
general depending upon the control objectives (mission objectives) and sensor/actuator 
locations. We will elaborate upon this in discussions below. 

It is therefore possible to produce reduced models in the form 

m r'% +k r<1r = e r u ’ (3.12) 

by employing the new method developed in Chapter n. In the representation of Eqn.(3.12), 
the variables qR represent the selected degrees of freedom that were determined to be of 
significance by user, and u is the input vector consisting of external control forces applied on 
the structure. 


3.2 Control Design 


Having presented a means of producing reduced models, it is now in order to 
investigate the role of these models in control design. Toward this task, we recall Eqn.(3.1) 
with the control term u introduced 

Mq+Kq = Eu ; qeR n , ueR m , (3.13a) 

and assume that the dof to be retained have been identified as q 2 and that the entire dof vector 
q has been reordered so that 

m * 

!< iil 


q - 


q 2 


q x eR ; q 2 eR 


(3.13b) 


We will assume that none of the dof contained in qi are directly acted upon by the control 
forces, so that the input matrix E has the form E T = [ 0 E 2 ]. With these considerations, 

Eqn.(3.12) in its partitioned form is written as 


' M 1 

M 12 

% 


K 1 

K 12 

V 


o' 

[ M n 

M2 

s 

+ 

[ K 12 

K 2 

q 2 


s 


u 


(3.13c) 


Then the reduced model, in terms of q 2 which are relabeled now as qR, is given by 

M R% + K R q R = e r u » (3.14a) 

where 
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(3.14b) 


MrAm^s,.* : Kr = V k L s i.2 • 

S u = (Kj-A/oc^-o ) 2 M 12 ). 


Theorem LL 

The reduced model (3.14) retains the frequency O) that was used to develop the model. 
Proof: Assuming u = 0 for the evaluation of natural frequencies, Eqn.(3. 14a) yields 

K “H| - | K 2 <2 S u-“ 2 < M 2- M n s .. 2 )| 

= |Kj - CO i’Mj - (Kj 2 - o^XK, - cb 2 M 1 ) ‘(K 12 - <o 2 M 12 ) | 


where Eqn.(3.14b), and the assumption that [Ki - oPMi] is nonsingular, have been used. 
Observing that the second determinant of Eqn.(3.15) yields the natural frequency of the 
structure concludes the proof. # 


1 


Kj - cu 2 ^ 


Kj - co 2 M 1 K 12 - to 2 ^ 

<2 ' Kj-CoX 


= 0 (3.15) 


An Illustrative Example. 


To illustrate the use of the reduced model (3.12) in control design, we use the simply 
supported beam shown in Fig.3.1 of Appendix A.. It is desired to control the linear 
displacement (W 2 ) and its linear velocity at node 2, by using the linear force actuator u located 
at node 2. The available sensors measure the linear position and linear velocity at node 2. An 
output feedback control policy is desired. 

For convenience of illustration, the beam is divided into two elements. Due to the 
statement of the problem, the W 2 dof is determined to be retained in the reduced model. Since 
W 2 is an intermediate dof, Eqn.(2.32c) from Variation B of the Algorithm, is used to develop 
the model : 


[s !j + s u ] l 


w„ 


0. 



(3.16) 
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from which 02 can be eliminated to yield 

S R w 2 = u . (3.17) 

By using ISrI = 0, the natural frequencies of the beam are extracted. These are {0.48, 2.15, 
5.32, 10.2 rad/sec.}. Once the natural frequencies are determined, a reduced model can be 
extracted from Eqn.(3.16) as 

M R (co) w 2 + K r (cd) w 2 = u, (3.18) 

where the argument CD is introduced to indicate that the corresponding parameters are 
dependent on the frequency at which Sr is computed; since there are four frequencies in this 
example, four reduced models can be generated. 

The model (3.18) is now used for control design. Though several control design 
schemes can be employed at this step, we have used the Linear Quadratic Regulator (LQR) 
theory [23]. Since the LQR theory requires a state-space representation of the model, 
Eqn.(3.18) is converted to 



Then the application of LQR theory yields the feedback control 

u = - g T x = - gj w 2 - gj w 2 (3.20) 

that minimizes 

oo 

J = J { w 2 (t) + W 2 (t) + p u 2 (t) } dt , (3.21) 

o 

subject to (3.19), where p is a design parameter, variation of which produces different gains 
{ g l » g2 } and hence, controllers of different bandwidth. Note that the control policy of 
Eqn.(3.20) can be implemented as an output feedback on the beam, since W 2 and W 2 are 
available as sensor signals. 

Note that this design methodology does yield an output feedback control policy as 
required, and it has been developed by min imizin g the objective function J of Eqn.(3.21) 
which represents the suppression of position and velocity at node 2. However, this policy 
has been designed based upon the approximation (3.18), and not based upon the true 
representation of the beam. Hence, in order to evaluate the performance of this control, it is 
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compared with the performance of the control that is designed based upon the true 
representation of the beam - a fourth order finite element model involving all the dof { 0 i, W 2 , 
02 , © 3 } has been used as a true model. The comparison is shown in Fig.3.2 of Appendix A. 

In Fig.3.2 of Appendix A, the labels 'regulation cost' and 'control cost' are defined 
as follows: 

00 

control cost = J u 2 (t) dt 

° - ( 3 . 22 ) 

A f 2 *2 

regulation cost = I { w 2 (t) + w 2 (t) } dt , 

0 

for some assumed initial conditions. Each value of p in (3.21) produces a point; each curve is 
produced by varying p and by calculating the control cost and regulation cost. The 'optimal' 
curve in this figure represents the performance of the system while using the control designed 
based upon the true model; the 'performance curve' corresponding to those controls that were 
developed based upon the approximation (3.18) would therefore be above the optimal curve. 
Since each natural frequency yields a different approximation (3.18), there are four 
performance curves in Fig.3.2 of Appendix A, each corresponding to one natural frequency. 

Observe that the approximation (3.18) obtained with © = ©1 yields the best 
suboptimal controllers, indicating that the lowest natural frequency is the most significant 
mode in this control problem. This is to be expected since the first mode has the maxim um 
effect on W 2 (the dof we wish to control) and on the sensors, and is affected the most by the 
actuator. Similarly the © = ©4 curve indicates that the highest mode is the least significant. 
Notice that the performance difference between the optimal curve and the © = ©1 curve is 
only marginal. The controllers from the controllers from © = ©1 curve are implementable as 
output feedback, while the optimal controllers (those from the optimal curve) are not 
implementable since they require feedback of the positions and velocities of all the dof, which 
in this example are not directly available. Therefore with © = © 1 , the approximation (3.18) 
produces fairly acceptable feedback designs. 

Another observation from this figure is that the closed loop system employing the 
suboptimal controllers remains asymptotically stable for all control cost. This is normally not 
expected of suboptimal controllers since the effect of errors in the approximation (3.18) is to 
drive the closed loop system unstable while using large feedback gains; see Ref.[24] for 
discussions on effect of model errors. This unexpected observation is explained by the 
following theorem. 
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Th eorem 32 

For a single input system, independent of the order of the full model (3.12), the feedback 
control designed employing the LQR theory based upon the approximation (3.14) is 
guranteed to yield the closed loop system atleast marginally stable, if the apprximation is a 
scalar equation (i.e., n 2 =l). 

The proof relies upon the following Lemma. 

LmmJLL 

The gains gi and g 2 involved in the control law 


U = -gl X -S2* 

(3.23a) 

that minimizes 


| { q t x 2 + q 2 x 2 + pu 2 } dt 
0 

(3.23b) 

subject to 


mx+ kx = eu , 

(3.23c) 

where m> 0,k>0, qj, q 2 , and p are arbitrary positive scalars, satisfy 

egj >0, eg 2 > 0. 

(3.24) 


Proof : 

Convert the minimizati on problem into the standard linear quadratic regulator problem as 


mm 


f 4 qi °lf x l 

f { [ X X ] T I + pu 2 } dt 

0 [° ^Jx 


subject to 


0 


± 0 
m 


[■] 


0 

e 

m 


u 


The solution to this problem is the full state feedback control [23] 


l r „ e, 
u = — 1° 

P m 


k k 
*11 *12 


12 


^22 


= “[^12 15**] 


where the ky (ij = 1,2) satisfy the Riccati equation 


(3.25a) 


(3.25b) 


(3.25c) 
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Comparing Equations (3.25c) and (3.23a) reveals that 
2 2 

e Sl=“ k 12 ’ e S2 = lk 22 * 

P P 


(3.25d) 


(3.26) 


Therefore, to prove Lemma 3.1, we need to show that k 22 > 0, and ki 2 > 0. 

Now, since qi and q 2 are positive, it follows from the positive definite property of the 
Riccati solution [23] that 

k^ > 0, k u > 0, k n k2 2 -kj 2 > 0. (3.27) 


Hecne, eg 2 > 0, as claimed in the lemma. To show that ki 2 > 0, rewrite Eqn.(3.25d) in 
its partitioned form as 


_k, 1 ,2 

- 2 ^ k 1 2 + < ll-— 2*12 
pm 

= o. 

(3.28a) 

k ll " 7^22 2^12^22 

pm 

= 0, 

(3.28b) 

2k .2- i ^4 + <l2 

pm 

= 0. 

(3.28c) 

Solve for ki 2 from Eqn.(3.28a) to get 



2k i2 , . 1,2 l' . 

= { -k+ / k + — } or 

pm V p 

{ -k- /k 2 +i ). 

V p 

(3.29a) 

Substituting the above result in Eqn.(3.28b) yields 


. ^ rr r 

k » = -sv - or 


(3.29b) 


Now, since kn must be positive, the correct solution for ki 2 is 
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(3.29c) 



Hence, egi >0. # 

Proof of Theorem 3.2: 

Since the approximation (3.14) is a scalar equation, it follows from Lemma 3.1 that the 
feedback gains in the resulting control policy 


u = *Si%-g24 R (3.30) 

satisfy E R gi > 0, and E R g 2 > 0. 

Now, since q R = q 2 and E R = E 2 , the control law (3.30), when implemented in the 
overall system of Eqn.(3. 13c) renders the closed loop system as 

M c q + D c q + K c q = 0, (3.31a) 


where 


M = M = 

C 


M 1 M 12 


M 


12 


M2 


> 0 



0 

0 



■ Kj K 12 

A* 

0 1 

E 2«2. 

K - 

’ *c 

K 12 K 2 +E 2 g l 

Now, since E 2 g 2 and E 2 gi are positive scalars, and since 



0 

0 



K = 

K + 


E 2 g l 



c 

0 




(3.3 lb\ 


(3.31b) 


(3.32) 


we get, D c > 0, Kc>0 (since K>0). Therefore, from the properties of mechanical 
systems[25], the system of Eqn.(3.31) is atleast marginally stable. # 

In conclusion of this Section, we point out that only a preliminary study of the 
applicability of the proposed method to control design has been presented. A more detailed 
study is in progress. The example included herein has served to illustrate the output 
feedback design using the new method. The results are satisfactory. An issue raised by this 
example is the following: Since each natural frequency yields a n approximation, one could 
conceivably obtain many approximations; as many as the total number of frequencies of the 
system. Hence, the development of the performance plot (as in Fig. 2 of Appendix A) 
including all of these approximations would be time consuming. This burden may be 
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reduced if the "significant" frequencies (such as the first frequency in the example) can be 
determined a priori - on emight adopt the modal cost anlysis of Ref. [5] for this task. 

The application of this method to a decentralized control design is given in 
Appendix B. 
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IV. CONCLUSIONS 

AND RECOMMENDATIONS FOR FUTURE RESEARCH 


The objective of this research effort was to develop a methodology that unifies the 
modeling problem and the control design problem. There exist a few methods in structural 
analysis and control design that may be integrated to produce such a unifying scheme. We 
have investigated the adaptability of the finite-element transfer-matrix (FETM) approach for 
this task. Specifically, we have developed a new method to overcome the limitations of the 
FETM method. This new method is capable of producing reduced models of structures, 
that can be tailored to suit the control design phase. The applicability of this method to 
controller design - both output feedback and decentralized control - has been illustrated 
through an example. 

The progress of this research has raised several issues that need to be addressed in 
future: 

1 . The new method requires an efficient numerical scheme for updating the assumed 
frequends during the extraction of the structural natural frequencies. An appropriate 
scheme needs to be selected or developed. 

2. The approximations produced by the new method, though yielding acceptabl 
controllers, do not preserve the symmetric properties of the original structure. The 
role of these approximations as reduced models is questionable. An appropriate 
modification of this method is hence required. This modeification is near 
completion and requires further investigation. 

3. Since the reduced model depends upon the frequency using which it has been 
derived, each natural frequency yields a different reduced model, and one could 
conceivably obtain many reduced models. This raises the question of which one (or 
which few) of the natural frequencies are to be used for the generation of the 
reduced models. In order to answer this question one needs to determine a priori the 
'significant' frequencies, such as those identified by modal cost analysis, for 
example. This issue has to be resolved for the proposed method to be a feasible 
controller design scheme. 

4. The controller design based upon the proposed method needs to be explored 
further. We believe that, since the reduced model is directly related to the full 
model, in terms of frequency, and the mass and stiffness matrices, the performance 
of the controller with respect to the full model could be predicted. In particular, the 
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stability of the overall closed loop system could be predicted, as established by 
Theorem 3.2 in a special case. The stability properties of matrix second order 
systems would naturally be used here. 

5. An important requirement of any control design is the robustness of the controller. 
This issue needs to be addresses in the future. 
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APPENDIX A. 


Control of Flexible Structures by 
Finite Element-Transfer Matrix Approach 

IF AC Symposium on Control of Distributed Parameter Systems, 
Los Angeles, CA, June 30 - July 03, 1986. 
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